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SUMMARY 


Power spectral densities of normal load factor have been obtained 
for two service operational training flights of a Republic F-8kG air- 
plane and three service operational training flights of a North American 
F-86A airplane in order to indicate the load-factor frequency content 
and possible uses of power spectral methods in analyzing maneuver load 
data. 


It was determined that the maneuvering load-factor time histories 
appeared to be described by a truncated normal distribution. 

The power spectral densities obtained were relatively level at fre- 
quencies below 0.05 cycle per second and varied inversely with approxi- 
mately the cube of the frequency at the higher frequencies. In general, 
the frequency content was very low above 0.2 cycle per second. 

The load-factor peak distributions were estimated fairly well from 
the spectrum analysis. In addition, peak load data obtained during 
service operations of fighter-type airplanes with flight time totaling 
about 2k, 000 hours were examined and appeared to agree reasonably well 
with the type of equations obtained from spectrum peak-load distributions. 


INTRODUCTION 


With a view toward providing information for the appraisal of the 
design requirements of future airplanes, the National Advisory Committee 
for Aeronautics with the cooperation of the U. S. Air Force and the 
Bureau of Aeronautics, Department of the Navy, has been conducting a 
statistical study of the pilot input, the airplane motions, and the 

-^Supersedes NACA Research Memorandum L56J15 by John P. Mayer 
and Harold A. Hamer, 1957- 
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loads imposed on a number of fighter-type airplanes during military 
service operations. The data obtained thus far are summarized in 
reference 1. 

The statistical treatment of the maneuver loads has been limited 
to the peak values of the loads and the frequency of occurrence. In 
the analysis of other loads associated with ;uch random factors as gusts, 
buffeting, or runway roughness (for example, refs. 2 to 5), the applica- 
tion of power spectral analysis methods has proved valuable, especially 
since a number of useful theoretical relationships such as the input- 
output relations and the frequency of occurrence of peak, loads have been 
established for stationary random processes. Such methods of analysis, 
however, have not been applied to maneuver loads. A single maneuver is 
clearly not a random process since it is associated with pilot intent. 
Nevertheless, maneuvering of a general nature (other than specific 
missions) over a long period of time or a collection of a great many 
maneuvers may acquire some of the characteristics of a random process. 
Therefore, it was believed that some of the relationships from the theo- 
ries of spectral analysis might be used as a guide in analyzing the 
maneuver load data. 

With this in mind, the normal load factors obtained in several 
flights of two jet fighter airplanes flown in regular squadron operations 
were analyzed by using power spectral methods. The data analyzed were 
obtained during the investigations reported : n reference 1 and they are 
a portion of the results summarized in reference 1. 

It is the purpose of this paper to prese nt the power spectra obtained 
in order to show the frequency content of normal load factor in service 
operations and to determine whether other ini ormation obtainable from 
spectral analysis methods (such as the frequency of occurrence of peak 
loads) could be obtained from the maneuver lead spectrum. 

Power spectral densities at frequencies up to 1 cycle per second 
are presented for two operational flights of a Republic F-84G airplane 
and for three operational flights of a North American F-86A airplane. 

The peak loads derived from these spectra are compared with the act ual 
peak loads. In addition, peak-load data frori other service airplanes 
are examined and compared with the results of spectrum analysis. 


SYMBOLS 


a 


f 


constant used in spectrum equation 
frequency, cps 
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•B 


"break" frequency, cps 

frequency of occurrence of normal load factor 

average number of normal load factor crossings per second 
of the zero axis with positive slope. 


f o 2 - 


I A(f) to df 

(ref. 1+) 


*( f W 


total number of positive normal- load-factor peaks per second, 

r ^( f w 


f 2 „ _±2 (ref. 4) 

P 


df 


m total number of observations 

N total number of positive peaks exceeding a given value 

n normal load factor, g units 

n(t) time history of normal load factor 

service limit normal load factor 

An incremental normal load factor (n - 1) 

Arij^ incremental service limit normal load factor ^n^ - lj 

Si mean incremental normal load factor 

An" filtered value of incremental normal load factor; at time t, 

Ant" = 0-9£nt_ At - 2Ant + 0.9Ant +At 

An q value of An at point of truncatxon 

p(y) probability that a given value y will be exceeded 
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P N (y) 

p i 

T 

T M 

t 

At 

a-(U 

&(f-0) 

S 


a 


®(f) 


An 


®< f >. An" 
X 2 


.05 


probability that a given value y will be exceeded based on 
a normal probability distribution 

probability of exceeding a given load-factor peak for load 
factors greater than 1 

total flight time 

maneuvering flight time 


time, sec 

time reading interval, sec 

proportion of area of normal curve represented by data in a 
truncated distribution 

Dirac delta function at f = 0 

mean of parent distribution from waich truncated distribution 
is derived 


root -mean- square value of An, a 2 = 


s(An - An) 2 


m 


'0 




power spectral density of incremental normal load factor (An), 

12 


g2/cps, «(f)^ n = 


lim i 
T— >» T 


r*T 

/ n 


j 


(t)i 


-2«ift 


dt 


-T 


(ref. 2) 


power spectral density of An” 
statistical reliability parameter 
critical statistical reliability parameter 


AIRPLANES 


The airplanes for which data were obtained were service models of 
the Republic F-84G and the North American F- }6 a. Both are low-wing 


L-1557 
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jet-propelled, fighter-type airplanes, the F- 86 A having a swept wing and 
empennage. The F- 86 A was equipped with an adjustable stabilizer and a 
hydraulically boosted elevator. 

The power spectra presented herein were obtained from flights made 
with external fuel tanks for the F-84G and without external fuel tanks 
for the F- 86 A. Except for the addition of sideslip and angle -of -attack 
booms, neither the external appearance nor the weight and balance of the 
airplanes was altered by the addition of the NACA instrumentation. 

Three -view drawings of the airplanes are presented in figure 1. 


INSTRUMENTATION AND TESTS 


The data were measured by standard NACA photographically recording 
instruments. In order to relieve the pilot of any recording-instrument 
switching procedure and thus to assist in obtaining normal operation, 
a pressure switch was employed to operate the recording instruments auto- 
matically at take-off. 

Normal load factors (normal to the longitudinal reference) were 
measured by an NACA air-damped recording accelerometer. The measuring 
element was damped to about O .65 of critical damping and the natural 
frequency was about 20 cycles per second. The accelerometer was located 
near the airplane center of gravity in such a way that the effects of 
angular velocities and angular accelerations were negligible. 

A standard two-cell pressure recorder connected to the airplane 
service system was used to measure the pressure altitude and indicated 
airspeed. The service systems were of the usual total-pressure -tube and 
flush static-pressure-orifice type. No corrections were made to the 
measurements for position error associated with the location of the 
static orifices. 

The flights were performed by service pilots during regular squadron 
operational training. Data were recorded continuously throughout a 
flight and were recorded only during those flights in which the mission 
was scheduled to include a large number of maneuvers. Although not 
requested, most of the maneuvers were performed in relatively smooth air. 

No attempt was made to specify the type or severity of maneuvers. Although 
the pilots were aware of the instrumentation, it was stressed that this 
was not to restrict their normal handling of the airplane since they would 
not be personally identified with the test results. 

The flights for which power spectra were obtained were selected since 
they contained most of the tactical maneuvers that are within the capa- 
bilities of the individual airplanes. A summary of the operations and 
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flight conditions for these flights is preset fted in the following table 
and includes the actual amounts of time from which the spectra were 
obtained. Distributions of the percentage of time (used for spectra) 
spent in various altitude and airspeed ranges for two combined flights 
of the F-84G and three combined flights of the F-86A are shown in fig- 
ures 2 and 3* respectively. 


Airplane 

Flight 

Pilot 

Mission 

Average 

weight, 

lb 

Avera^ e 
pressi re 
altituc e, 
ft 

Average 

indicated 

airspeed, 

knots 

Average 

Mach 

number 

Time for which 
spectrum was 
obtained, sec 

Republic 

F- 81 +G 

1 

2 

A 

B 

Acrobatics and 
dive-bombing 

16,000 

11,00( 

290 

0.33 

3,900 

Acrobatics, 
dive-bombing, 
and ground 
gunnery 

2,777 

North American 
F-86A 

1 

C 

Acrobatics 

13,060 

1 A, 15 ( 

298 

00 

1,876 

2 

C 

Acrobatics 

1,823 

3 

1 D 

Dogfighting 

1,032 


METHODS OF ANALYSIS 


Power Spectral Densities 

In the present paper use will be made of some of the concepts from 
the theories and methods of power spectral analysis. The methods used 
in the present paper are described in referex.ces 2 and 3* 

The methods described in reference 2 were used in determining the 
power -spectral densities of the normal -load -factor data. The records 
were read at 0.5-second intervals for each flight and 40 estimates of 
the power at frequencies up to 1 cycle per second were obtained by using 
digital computing methods. In order to improve the accuracy of the com- 
puted spectra, the spectra were first obtained of an altered (filtered) 
time history of load factor given by the equation 


Ant" = °*9 ^t-At - 2 An t + 0.9 ^t +A t 


( 1 ) 
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This is essentially a high-pass filter which retains some information 
at zero frequency. This particular expression was used since it was 
found to make the spectrum of An" relatively constant over the fre- 
quency range . The spectrum of An" was then converted to the spectrum 
of An by the equation 


4>(f) 


®(f) 


An' 


^ (1.8 cos 2*f At - 2) 2 


(2) 


In addition to the spectra computed from the records read at 
0 . 5-second intervals, spectra were also computed from the data taken 
at 2-second intervals and at 10-second intervals in order to provide 
more detail on estimates of power at the lower frequencies. The spec- 
tra computed for 2-second intervals were computed for each flight at 
frequencies up to l/4 cycle per second. The spectra computed for 
10-second intervals were computed for the combined individual flights 
for each airplane at frequencies up to l/20 cycle per second. The data 
at 10-second intervals were not filtered by using equation (l) since 
the spectrum was relatively level. 


Ordinarily spectra are obtained of quantities which fluctuate about 
some mean value. Maneuver load time histories, on the other hand, fluc- 
tuate about 1 g but generally only in the positive direction. 

The power spectral densities of the incremental maneuvering load 
with a nonzero mean (An) would include a Dirac delta function at the 

origin with an area (£n) : 

oCf)^ = 2A^ 2 6(f-0) + ^(f)^ (3) 

The power spectra presented in this paper have the mean taken 
out and are equal to the last term oCf)^ in equation (3). 


Load-Factor Counts 

Threshold counts . - The probability curves for the normal-load- 
factor data were determined by a process which minimized the level- 
flight portions of the record. In this method counts were made at a 
given threshold every time the record crossed the given threshold and 
exceeded one -half of the increment to the next threshold. In this paper 
threshold values of 0-5g were used. For example, counts were made at 
2g if the load-factor record crossed 2g while increasing and then 
increased beyond 2.25g or if the record crossed 2g while decreasing 
and decreased beyond 1.756* This method is illustrated in figure 4 
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(method 5)* In using this method all small fluctuations in the load 
factor less than + 0 . 25 g are eliminated. 


Peak counts . - Two methods of obtaining major peak distributions 
were used in this paper and are illustrated in figure 4. In the first 
method (method l) peaks were counted with a threshold of 0.25An or 
0 . 25 g whichever was greater; that is, the load factor must decrease 
more than one— quarter of the value of the peak— load— factor increment 
An before and after the peak. For values of An below 1 (2g), how- 
ever, a threshold of 0 . 25 g was used in which the load factor had. to 
decrease 0.25g before and after the peak. The second method (method 2) 
was method B of reference 1. In both methods peaks of less than 0.25g 
were not counted. Peaks obtained by using mei hod 2 were counted manu- 
ally. Those obtained by using method 1 were counted automatically from 
the time histories with a digital computer. 


RESULTS AND DISCUSSION 


Spectrum for F-84G Airplane 

In figure 5 the power spectrum of normal load factor is shown for 
flight 1 for the F-84G airplane. Results are shown for the four quar- 
ters of the flight, of which each part consisted of about 900 seconds 
of flight. The root -mean-square value a and the mean 251 for each 
part are shown in the figure. 

The power spectral densities for each of the four parts appear to 
be very similar and decrease approximately with f-3 . 

The total spectrum for flight 1 of the F-34G airplane is shown in 
figure 6. Also included is the spectrum of the data taken at 2-second 
intervals to provide more detail at the lower frequencies. In addition, 
the statistical 95-percent confidence bands associated with the spectrum 
are shown. For simplicity, the confidence ban is shown were obtained 
for a faired curve rather than for each indivilual point. The value of 
°( f )An at zero frequency shown has little statistical reliability 
and is shown only as a matter of interest. 

The confidence bands shown in figure 6 are representative of those 
obtained for all the spectra presented in this paper. Therefore, con- 
fidence bands for the other spectra will not be given. 

The spectrum obtained with At = 2 second;; provides more detail 
at the lowest frequencies and generally agrees with the spectrum obtained 
with At = 1/2 second except at the highest frequencies . At the highest 
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frequencies for the 2 -second data (i.e., l /4 cps) the spectrum is not 
expected to be accurate since the power at higher frequencies is included 
in the spectrum obtained. Therefore the values for the 2-second data 
are not shown at the higher frequencies. 

Effect of level flight on spectrum . - The effect of the presence in 
the time history of long intervals at 1 g or level flight on the spec- 
trum is shown in figure 7 . The results shown are for about 900 seconds 
from flight 2. The record was examined and approximately 530 seconds 
of flight near 1 g were deleted. The resulting record was then analyzed 
with one maneuver connected to the end of the previous one. The results 
are presented as the power spectral densities divided by the mean square 
CT 2 of the respective parts. The value of 0 2 of the original time his- 
tory was 0.398 and the value of a 2 for the all -maneuvering time history 
(original less level flight) was O. 929 . It may be seen from figure 7 that 
above 0.025 cycle per second the spectra are directly proportional to the 
mean square or 


where 

flight 


% r. 
time. 






L ° 2 

1 T 

a 2 

Tm 

maneuvering 

h 

flight 

time 

o 2 t 

% _ 

> (f U T 


T 



T 


( 4 ) 


to the total 


( 5 ) 


In general, the effect of level flight on the spectrum is to change the 
level of the spectrum. 

Spectra for flights 1 and 2 .- In figure 8 the spectra for the com- 
plete flights 1 and 2 are given. In each case the power has been divided 
by the mean square for the particular flight. The power spectral den- 
sities are seen to be quite similar for the two flights. At the higher 
frequencies (above 0.1 cps) the two curves have the same shape but are 
of different magnitudes. 

In figure 9* flights 1 and 2 have been combined and, in addition, 
the power spectral densities obtained for the combined flights 1 and 2 
at 10-second intervals have been included. The 10-second data have been 
corrected for "folding" based on the spectra obtained from the l/ 2 -second 
and 2- second data. When the 10- second data are included, a peak appears 
at about 0,025 cycle per second which corresponds to the average number 
of maneuvers per second during the two flights, which is roughly 0.024 
(86 maneuvers per hour). 
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Spectra for F-86A Airplane 

The power spectral densities for three operational flights of the 
F-86A airplane are shown in figure 10. The results are ratioed to the 
mean square values of the incremental load fac <or . The spectra for 
these three operational flights also appear to he similar. 

The power spectral densities for the comb .ned flights 1, 2, and 3 
for the F-86 airplane are shown in figure 11. Included in this figure 
are the data obtained at 10-second intervals. As was the case for the 
F-84G airplane, there appears to be an indication of a peak in the spec- 
trum at approximately 0.02 cycle per second. Again this value corre- 
sponds to the average number of maneuvers per second in flights 1, 2, 
and 3 of about 0.022 (79 maneuvers per hour) . 


Comparison of Spectra for F-84G and F-86A Airplanes 

In figure 12 are shown the normalized power spectral densities for 
three combined flights for the F-86A airplane and two combined flights 
of the F-84G airplane. The same curves are shown on a linear scale in 
figure 13 . 

The power spectral densities for the two airplanes in squadron 
operations are very much alike (fig. 12); however, there appears to be 
some increase In power at the highest frequencies (l cps) for the 
F-86 airplane as compared with the F-84G airplane. It is believed that 
this increase is caused by a lightly damped longitudinal oscillation 
which was present in the F-86 airplane. The effect of this is probably 
present at frequencies between 0.4 and 1 cycle per second. In figure 12 
the power spectral densities are shown to be relatively constant at fre- 
quencies below 0.03 cycle per second. The pover is very low above 
0.2 cycle per second. (See fig. 13*) 


Analytical Representation of Maneuver Load Spectrum 


The spectra of figure 12 appear to be of the type given by the 
formula for the simple spectrum: 


«(fW 


1 + ( f / f B ) a 


(6) 


where fg is the ,f break" frequency. In figure l4 the total spectra for 
the F-84G and F-86A airplanes (shown in fig. 12) were averaged according 
to the flight time for each airplane and are compared with equation (6) . 
In this case a was found to be 2.6 and f B was 0.031 cycle per second. 
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This equation appears to fit the data above a frequency of 0.005 cycle 
per second. Also shown in the figure as a matter of interest is the 
approximate phugoid frequency range which is between 0.008 and 0.025 cycle 
per second. As previously noted, the average number of maneuvers per sec- 
ond was 0.022 and 0.024 per second for these airplanes (79-86 per hour). 


Probability Distributions 

Threshold distributions .- The frequency distribution of threshold 
counts of normal load factor for the combined three flights of the 
F-86A airplane is shown in figure 15 and table I. In examining the 
frequency distributions, it was noted that the distributions looked 
much like normal distributions if only the positive or only the nega- 
tive distributions were examined. For example, in figure 15, normal- 
distribution curves are shown for both the positive and the negative 
distributions. These curves may be described as truncated normal dis- 
tributions above and below An = 0. (See ref. 6 for a discussion of 
truncated distributions.) Since the distributions are essentially one 
sided, the measured value of the threshold counts at An = 0 must be 
adjusted. Therefore, in figure 15 it may be noted that the zero value 
for the positive side is given as two times the number of zero threshold 
values having a positive slope, and. In a similar manner, the zero value 
for the negative side Is given as two times the number of zero threshold 
values having a negative slope. In table I the threshold counts labeled 0 
are the zero values having positive slope and those labeled (-)0 are the 
zero values having negative slope. 

The probability of exceeding a given incremental normal load 
factor An^ r for a truncated normal distribution is 



for Antr ^ 0 
where 


Ant r = An - Ano 


and 

A^r - 5 


Z 


a 



1J 


The normal -load -factor data are truncated at An = 0 and the mean 
of the parent distribution is assumed to be at An = 0. Therefore, 

= 0.50 

= 0 


l/Wi \ 2 

e df£\ An > 0 (8) 

\ a J 

Since it is sometimes convenient to examine data on normal proba- 
bility graph paper, the quantity a( ^)P(An) of a truncated normal dis- 
tribution can be plotted and the results will fall in a straight line 
on normal probability paper. 

The probability distributions of threshold counts of normal load 
factor for flights 1 and 2 of the F-84G airplane and flights 1, 2, and 3 
of the F-86A airplane are shown on a normal probability scale in fig- 
ures lb and 17 . The data are shown for the individual flights as well 
as for the combined flights for each airplane. The value a(f;)P(An) 
is plotted against An. In each case, the zero value used was two times 
the number of zero thresholds having a positive slope. The threshold 
counts for each airplane are given in table I. The points shown in 
figures lb and 17 are plotted at l/4g less than the threshold values 
used (i.e., 1 plotted at 0.75, 1.5 plotted at -. 25 , and so forth) since 
the thresholds were based on bands ±l/4g about the threshold values. 

In figures lb and 17 it may be seen that nhe truncated normal dis- 
tributions appear to fit the data. In addition to the plots on the 
normal probability scale the (Chi-square) statistical test was made 
for each flight and the combined flights. (Se ; ref. 7, for example.) 

The results of these tests as shown in the figures are well within the 
95-percent level of significance and indicate that the maneuvering loads 
of the operational flights presented in this paper can be considered as 
truncated normal distributions. 

Although similar analysis could be made of the negative maneuvering 
load factors, the data sample sizes used were not sufficient for any 
reliable analysis to be made. 

Peak-load distribution . - The peak -load -fa :tor frequency distribu- 
tions obtained for the combined flights of each of the airplanes are 
given in table I. Two peak distributions are given: the total peak 

distribution and the major peak distribution. The total peak distribu- 
tion includes all peaks obtainable from the tine histories of normal 


a( 0 

i 

and 

O oo 

P(An) = 

0 • 5 y 2it j 

U An/cr 


-o \j] vji h tr -1 
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load factor. Since the load factors were calculated to O.OOlg, the 
threshold for the total distribution is O.OOlg, which is less than the 
reading accuracy of about O.Olg. A count was made of all the peaks by 
using a threshold value of O.Olg for the F-84G airplane and it was found 
that the difference between the frequencies counted by using the two 
thresholds was almost entirely included in the 0 to 0.25g range. There- 
fore, the values beyond 0 . 25 g are probably not affected to any extent 
by the threshold used. In addition, since the records were read at 
0. 5-second intervals the minor peaks at high frequencies were also fil- 
tered out; however, these are not believed to be of any importance in 
these flights. 

In figures 18 and 19 the major peak cumulative frequency distribu- 
tions are shown plotted against (An) 2 for the F-84G and F-86A airplanes. 
Also plotted in figures 18 and 19 are the threshold frequency distribu- 
tions. From equation (8) it may be seen that the threshold frequencies 
should plot as a straight line when plotted against (An) 2 on semilog 
paper if the distribution is a truncated normal distribution. Also it 
may be seen from figures 18 and 19 that the major peak cumulative dis- 
tribution plots approximately as a straight line against (An) 2 above 
An = 1. 


Calculated peak -load distributions .- In reference 4 an expression 
for the probability of exceeding a given peak is given for a Gaussian 
random process. Although the maneuvering load data of this paper can- 
not be considered a Gaussian random process, the results of figures l6 
and 17 indicate that the positive maneuver loads can be considered as a 
truncated or "half" normal distribution. Therefore, it was believed 
that perhaps the peak probability expression given in reference 4 might 
be used as a guide in calculating the peak distributions obtained in 
the present tests. 


If the maneuvering loads were Gaussian in character the probability 
that a peak load will exceed a given value can be given, as in refer- 
ence 4, as 


P(An) 



(9) 


where 

Kl = \/l - (fo/fp) 2 


and 


k 2 


Kl 

f o/% 
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and where the normal probability distribution Pjj(y) is 



Poo 

_ l 'An_\ 2 

/An \ 

_ 1 . / 

e 2 d/— 2— \ 

IkW 

\/2 « / An 

M 


Kpa 


( 10 ) 


where i = 1 or 2. 

The total number of positive peaks exceeding a given value is 

N = f p T M P(An) (11) 

For large values of Ai equation (9) approaches 

_ l/fe] 2 

P(An)*^e 2 ^ a) (12) 

f P 


and for the asymptotic case the number of pea* loads exceeding a given 
value approaches 


_i(An') 2 ( 13 ) 

N - f 0 T M e / 

Since much of the flight time is spent in nonmaneuvering or near- 
level flight, the maneuvering flight time rather than the total flight 
time is used in equations (ll) and ( 13 )* In order to determine the 
maneuvering flight time, the percentage of total flight time spent above 
a given load factor was determined. These times are shown in figure 20 
for the F-84G airplane and in figure 21 for the F-86A airplane. The 
distribution of the time spent in various load-factor intervals is given 
in table II. Since in the threshold counts and in the peak counts the 
minimum load-factor variation was 0 . 25 g, it was assumed that the maneu- 
vering flight time was the portion spent above An = 0.25g. It can be 
seen in figures 20 and 21 that about 52 percent of the total flight 
time could be considered maneuvering flight time for both the F-84G air- 
plane and the F-86A airplane in these particular flights. 

The peak -load data for the F-84G airplare are shown in figure 22 
and the peak-load data for the F-86A airplane are shown in figure 25 . 

Two peak-load distributions are shown - the total peak distribution 
and the major peak distribution. The major peaks shown in figures 22 
and 25 were counted by both method 1 and method 2 of figure 4; however, 
only the results from method 1 are shown. Tie results from method 2 
were almost the same as those from method 1. 
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The peak-load distributions as given by equations (ll) and ( 13 ) 
are also shown in figures 22 and 23- The value of the root -mean-square 
load factor o used was that obtained from the threshold analysis of 
figures lo and 17 . The value of fo/^p obtained from the integration 
of the spectrums was 0.184 for the F-84G airplane and 0.135 for the 
F- 86 A airplane. The value of fQ was O.O 76 for the F-84G airplane and 
0.08l for the F- 86 A airplane. 

It may be seen in figures 22 and 23 that both the total peak count 
and the major peak count are approximately estimated by equations (ll) 
and ( 13 ), respectively. The agreement seems somewhat better for the 
F- 86 A airplane than for the F-84G airplane. 


Although the total peak count is of interest, it is the major peak 
distribution which is of most value. For example, in fatigue analysis, 
it is believed that the major peak distribution would more accurately 
describe the life history of loads than would the total peak distributions . 


Major peak-load distribution in service operations . - In order to 
compare the results from the spectrum analysis with more complete flight 
data from regular service operations the major load-factor peaks obtained 
in the service operations of references 1 and 8 to 11 are shown in fig- 
ures 24 to 26 . The airplanes for which information is presented in 
figures 24 to 26 are 


Republic F-84E and G 

North American F- 86 A, E, and F 

North American T-28A 

North American B-45A 

North American AJ -1 

Lockheed F-80A, B, and C 

Lockheed F-94B 


Lockheed TV -1 
Lockheed P 2 V -3 
Lockheed P2V-4 
Grumman F 8 F -2 
Grumman F9F-2B 
McDonnell F2H-2 
Douglas AD -4 


All data in these figures were obtained during training operations unless 
otherwise stated. The data represent about 24,000 hours for f ighter-tvpe 
airplanes and about 2,000 hours for bomber-type airplanes. In order to 
prepare figure 24 the total number of load-factor peaks exceeding a given 
load factor N were plotted against the parameter (An/An ^) 2 f or 

airplane. A straight line was then fitted to the data above a value of 
(An/An L j 2 of about 0.2, and the straight line extrapolated to the value 
N = No' at zero-load-factor increment. The ratio N/n 0 * of the number 
of load-factor peaks exceeding a given load factor to the extrapolated 
value of exceeding a value of An/Anp =0, is plotted for each airplane 
in figure 24. The solid symbols in figure 24 represent less than 10 points 
and are considered less representative of the peak distributions than the 
other points. Also shown in figure 24 are the number of hours of flight 
represented by the data and the value of Nq'/T* the extrapolated value 
of the number of peaks to exceed 1 g per hour. 
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For large values of load factor the peak 
the form of equation ( 13 ): 





distribution is given in 


(14) 


where 


Nr 


f t Tm 
= lo 1_ T 


(15) 


and 


N 0 /T = fo^ (16) 

When plotted on semilog paper as in figure 24, equation (l4) is a 
straight line with the slope - - 4 -. The values of the root mean square a 

2 a 2 

for the data shown vary between 0,25 and 0.55- The average value of a 
in figure 24 is 0.284 and it appears that the straight line fits most of 
the data fairly well above lo ^£n/Anp^ ^ = 0 . 0806 ^ at least up to the 

limit load factor. Of course, the value of <i is to some extent a func- 
tion of the mission of the airplane and it is possible that, for a spe- 
cific mission, considerably different values than those shown in figure 24 
could be obtained. For example, in figure 24, the training data for the 
F-84E and G airplanes were obtained in Korean operations in which the 
airplane was used mostly in air-to-ground missions. The value of the 
root mean square for these data appears to be considerably lower than 
that for the other airplanes. (This effect would also be apparent if the 
true service limit load factor was lower than that assumed ( 7 - 33 ) in fig- 
use 24.) Most of the other training data shown in figure 24 for the 
other airplanes were obtained in training whi :h included many types of 
missions . 

It may be noted from equation (l 6 ) that, if it can be assumed that 
f 0 is about 0.08 and does not change, the percent of maneuvering flight 
time can be determined from the values of Nq/t given in figure 24. 

3d = nq'/t = V/T (l7) 

T f Q (0 .08) (5601) 

From the values of Nq'/T in figure 24 and from equation ( 17 ) it was 
determined that approximately 2 to 8 percent of the total flight time 
was spent in maneuvering in the Air Force training operations, about 
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2 to 6 percent was spent in maneuvering in combat operations in Korea, 
and between 7 and 24 percent in the NACA tests in operational training. 
The value of 52 -percent maneuvering -flight time obtained for the five 
flights of figures 20 and 21 are applicable to just those particular 
flights which were more active than the average. This percentage is 
obviously a function of the airplane mission and would be higher for 
gunnery or dive bombing training than for other operations. 

In addition to the data available from Air Force operations there 
are data obtained with V-G recorders from U. S. Navy operations. (See 
ref. 11.) These data are presented in figure 25 as plots of the time 
to exceed 1 g in maneuvers to' to the time to exceed a given load 
factor t. Since, in V-G records a true frequency count cannot be 
obtained at low load factors, the data were used only where the count 
appeared to be approaching the true frequency. The times to exceed a 
given load factor were plotted against ^An/ ^ an< ^ a straight line 
was fitted to the data and extrapolated to zero -load -factor increment 
and then the ratio of this extrapolated zero time to to the time to 
exceed a given load factor was plotted in figure 25- It may be seen 
that 


V _ _N_ 
t ‘ N q 1 


(18) 


The straight line shown in figure 25 is that obtained from figure 24 
with a value of o = 0.284. The black symbols represent points which 
consist of less than 10 peaks. The line appears to fit the data fairly 
well for most of the airplanes shown. 


In figure 26, data obtained from references 10 and 11 are shown 
for several large airplanes . Also shown is the line obtained from the 
data for the fighter-type airplanes of figures 24 and 25- It may be 
seen that the data for the AJ-1 and P2V-4 airplanes are close to the 
line; however, the data from the B-45 and P2V-3 airplanes appear to have 
a higher root -mean-square value than that of the line shown. It is evi- 
dent that the line shown would probably not be a good estimation of the 
loads to be encountered in airplanes which are primarily nonmaneuvering 
airplanes such as large bombers, but would be of more use in more maneu- 
verable airp lan es. In addition,, it is not believed that the results for 
other low-limit load-factor airplanes would necessarily indicate more 
severe use (as indicated for the B-45 and P2V-3 airplanes) than the line 
shown. It is believed that for any large bomber the type of mission 
would have a large effect on the slope of the line. 

Other methods of treating load- factor peak distributions .- From the 
spectrum analysis it was indicated that the major peak-load distribution 
could be given by equation ( 13 ) above la and this is also indicated 
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from figures 2k and 25 . It can be seen that below la, however, the 
peak -load distribution is not well represented by equation ( 13 ). One 
way of accounting for the peaks at low load factors is to assume that 
the major peak distribution is made up of two peak distributions having 
different values of a. If such is the case, the probability of exceeding 
a load -fact or peak greater than 1 g can be given as 


Pi = 


N 


f 0 T M 


3e 

t m 


/An/£n L \ 2 

l/An/An L \ 

l 01 J +% 

2 \ J 


M 


where 


T^_ maneuvering time associated with o- j_ 
T 2 maneuvering time associated with a 2 
For the data of figure 2k, 


( 19 ) 


o 1 = 0 . 281 + 


a 2 = 0 * 100 
T-^T = T 2 ^T = 0.50 

In figure 27 the data of figure 2l+ are shown with the probabilities 
given by equation ( 19 ). The probability given by equation ( 19 ) is 
almost identical to the "Standard Probability Curve for Maneuvers" given 
in reference 1. It may be seen that the data are represented fairly 
well by equation ( 19 ). 


In figures 22 and 23 it might be noted that, although the major 
peaks were approximated by equation ( 13 ), the curve from the equation 
which represented the total peak count appears to be similar in shape 
to the major peak count. Therefore the major peak count could be approxi- 
mated by an equation similar to equation ( 9 )- 


The number of load -factor peaks exceeding a given load factor 
divided by the number of load -factor peaks exceeding 1 g from equa- 
tion (9) is 


N 


N, 


'0 



1 + f o/ f P 


■N 


/An/An L ' 


K-jO 





1 



(20) 


L-1957 
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When n/Nq is plotted against y — L — Lj it is approximately linear 

above 2 a and has a zero intercept of about ip q (at least for the 
present case of f Q /f p = Therefore, if the total number of 

major peaks above a given load factor are plotted against ^Any^n^)^ 
and a straight line fitted to the data, the slope of the line will 
approximately determine the root-mean-square value a. Taking into 
account that the straight-line extrapolation of equation (20) will have 
a zero intercept of IPq or INq the experimental data such as those 

in figure 2k- can be compared with equation (20) if it is plotted as 

N/2Nq. In this form N/Nq = p/Pq = P 1 which is the probability of 
exceeding a given load factor for load factors greater than 1 g. In 
figure 27 the data are compared with equation (20) in this manner and 
it is in fairly good agreement. 

Thus, any of the three equations (13), (l9)> and (20) can be used 
to approximate the major peak -load probability curve. Probably equa- 
tion (13) is the simplest if the shape of the curve below la is unim- 
portant. If the low load factors are important, for instance in some 

fatigue applications, equation (19) or (20) would be more desirable. 

In addition, probably none of these equations are valid much above the 
limit load factor since inadvertent maneuvers may affect the curves 
mostly at high load factors. (See refs. 1 and 12.) 

CONCLUSIONS 


The results of the power spectral analysis of the maneuvering load 
data indicate that spectral methods appear to have some promise in the 
statistical treatment of maneuver loads data. It is admittedly an 
almost empirical use of these methods; however, it is an attempt to 
describe more completely the complex process of maneuvering loads on 
airplanes other than by simple peak counts. In addition, it should be 
pointed out that the flights for which the spectra were obtained were 
training flights of a general nature in which several different missions 
were simulated and that the results for an airplane flying one specific 
mission could be much different. From the present analysis the following 
conclusions have been reached. 

1. The maneuver-load-factor time histories of the five operational 
flights examined could be described by a truncated normal distribution. 

2. The power spectral densities obtained were relatively constant 
at frequencies below 0.03 cycle per second and varied inversely with 
approximately the cube of the frequency at the higher frequencies up to 
1 cycle per second. Above a frequency of 0.005 cycle per second the 
spectrum could be approximated by a simple spectrum equation. 
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3. The frequency content above 0.2 cycle per second (5-second 
period) was very low. 

4. The load -factor peak distributions were estimated fairly well 
from the spectral analysis. 

5. Peak-load data from about 24,000 hours of service operations 
of fighter -type airplanes appeared to agree reasonably well below the 
limit load factor with the type of equations obtained from spectrum 
analysis . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 2, 195 b * 


L-1557 
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TABLE I.- FREQUENCY DISTRIBUTION OF NORMAL LOAD FACTOR 
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TABLE II.- FREQUENCY DISTRIBUTION OF TIME SPENT IN 
VARIOUS NORMAL-LOAD-FACTOR INTERVALS 


Time, sec 
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Figure 5 ■ - Power spectral densities of normal load factor for four parts 

of flight 1, F-84G airplane. 
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Figure 6.- Power spectral densities of normal load factor for flight 1, 

F-84G airplane . 








Figure 9*- Power spectral densities of normal load factor for two coim 
bined operational flights of the F-84G airplane . 
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combined operational flights of the F-86A airplane . 




Figure 12.- Comparison of power spectral densities of normal load factor 
for two operational flights of the F-84G airplane with those for 
three operational flights of the F-86A airplane . 




1 

O F-34G 

, Flights 

1 and 2 

□ F-36A 

9 Flights 

l,2 # and 3 


lull 


Frequency, c ps 


Figure 13.- Comparison of power spectral densities of normal load factor 
on a linear scale for two operational flights of the F-84G airplane 
and three operational flights of the F-86A airplane. 
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Figure 15.- Illustration of truncated normal distribution used in 
describing threshold counts of normal load factor for F- 86 A air- 
plane, flights 1 , 2 , and 3 . 
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(a) Flight 1. 


(b) Flight 2. 



(c) Flight 3 - 


(d) flights 1, 2 , and 3« 


Figure 17.- Probability that a given threshofd value of load factor will 

be exceeded. F-86A airplane. 
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Figure 20.- Amount of flight time spent above given load factors for 
flights 1 and 2 , F-84G airplane. 
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Figure 21.- Amount of flight time spent above given load factors for 
flights 1, 2, and 3, F-86A airplane. 
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Figure 22.- Comparison of load factor peak distributions with those 
obtained from power spectral densities. F-84G airplane. 
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Figure 23-- Comparison of load factor peak distributions with those 
obtained from power spectral densities. 7-86A airplane. 
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Figure 24.- Percentage of load factor peaks exceeding a given value 
°f (An^/An/j obtained in squadron operations for fighter-type 
airplanes . 
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Figure 25.- Percentage of load factor peaks exceeding a given value 
of ^An/An^ obtained from V-g records obtained during squadron 
operations for fighter-type airplanes. 





Figure 26.- Percentage of load factor peaks exceeding a given value 
of ( An /zn T ^ ^ obtained in squadron operations for large airplanes 
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Figure 27.- Comparison of load factor peak cumulative distributions 
with analytical distributions for fighter-type airplanes. 


NASA - Langley Field, Va. Ti — 1 



